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ABSTRACT 

We introduce exact analytical solutions of the steady-state hydrodynamic 
equations of scale-free, self-gravitating gaseous disks with flat rotation curves. 
We express the velocity field in terms of a stream function and obtain a third- 
order ordinary differential equation (ODE) for the angular part of the stream 
function. We present the closed-form solutions of the obtained ODE and con- 
struct hydrodynamical counterparts of the power-law and elliptic disks, for which 
self-consistent stellar dynamical models are known. We show that the kinematics 
of the Large Magellanic Cloud can well be explained by our findings for scale-free 
elliptic disks. 

Subject headings: hydrodynamics — galaxies: kinematics and dynamics — galax- 
ies: structure — methods: analytical 
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1. INTRODUCTION 

Astrophysical disks are either stellar or gaseous. The existence of non-axisymmetric 
equilibrium states of such systems has been attractive both from theoretical and observa- 
tional aspects. The search for dynamical equilibria of stellar systems is known as the self- 
consistency problem where one seeks for a positive phase space distribution function (DF) 
that produces the density of the model and satisfies the coUisionless Boltzmann equation 
(Binney & Tremaine 1987). The construction of the DF can be done by analytical means or 
numerical exploration based on Schwarzschild's orbit superposition method (Schwarzschild 
1979; Kuijken 1993, hereafter K93; Levine & Sparke 1998; Teuben 1987; de Zeeuw, Hunter, 
&; Schwarzschild 1987; Jalali & de Zeeuw 2002, hereafter JZ). The main difficulty linked 
to the self-consistency problem is the lack of the second integral of motion for most planar 
models. 

Gaseous disks, however, are investigated through solving the hydrodynamic equations of 
compressible, self-gravitating fluids. There are several works in the literature that deal with 
the stability of perturbed axisymmetric gaseous disks (Heemskerk, Papaloizou, & Savonije 
1992; Goodman & Evans 1999; Lemos, Kalnajs, & Lynden-Bell 1991). The problem of non- 
axisymmetric systems is somewhat different because our knowledge is limited even about 
equilibrium states. Following a first-order hnear theory, Syer & Tremaine (1996, hereafter 
ST96) constructed simple scale-free disks and approximately solved the steady-state hy- 
drodynamic equations for the components of the velocity vector. They also proposed a 
nonlinear theory based on the calculation of Fourier coefficients using Newton's method. A 
similar strategy was adopted by Galli et al. (2001, hereafter GOl) in the study of isother- 
mal molecular clouds and the fragmentation hypothesis of multiple stars. Their nonlinear 
solutions, in contrast with ST96, were generated by direct numerical integration of reduced 
ordinary differential equations (ODE) and an iterative scheme to match the boundary con- 
ditions. These methods of approach seem inefficient when one seeks for global properties of 
highly non-axisymmetric equilibria and the transient response of fluid disks in the vicinity 
of such states. Therefore, flnding analytical solutions for the hydrodynamic equations is of 
special interest even in the steady-state conditions. 

In this paper we explore exact non-axisymmetric solutions of self-similar gaseous disks 
with flat rotation curves. In §2 we express the governing equations and reduce them to 
a third-order ODE in §3 where we also introduce our closed-form analytic solutions. We 
construct hydrodynamical counterparts of the JZ models in §4 and apply our results to the 
Large Magellanic Cloud (LMC) in §5. We show that the kinematics of the LMC can be 
explained by our model elliptic disks. We end up with conclusions in §6. 



-3- 



2. STEADY-STATE HYDRODYNAMIC EQUATIONS 

We consider a thin disk of compressible, inviscid self-gravitating fluid in a steady 
state (dynamical equilibrium). The hydrodynamic equations can be written as the time- 
independent continuity equation 



and the momentum equations 



^9 1 9 2N , ^vrv^ 1 dp 1 dU 

-^-Q^mvnv,) + -^(E.,) + ^ = --_ - E--, (3) 

where {R, 0) are the polar coordinates, v = {vn, v^) is the velocity field in Eulerian descrip- 
tion, and p is the pressure. The potential U is due to the self-gravity, which is related to the 
surface density E through the Poisson integral 

U(R6)- C ff m'A')R'dR'd<j,' 

J J y/R^ + R'^ - 2RR' cos{4) - 0') ' 

where G is the gravitational constant. In the following section, we attempt to simplify the 
governing equations for scale-free disks and solve the resulting equations by analytical means. 



3. NON-AXISYMMETRIC DISKS WITH LOGARITHMIC POTENTIALS 

Gaseous disks with logarithmic potentials possess the astrophysically important prop- 
erty of a flat rotation curve. The circular gaseous disks having this property are known as 
Mestel (1963) disks, the stability of which was investigated by Goodman & Evans (1999). 
Here we derive exact analytical expressions for the physical quantities of non-axisymmetric 
disks. 



3.1. Stream Function 



We assume cuspy, scale-free potential-density pairs of the form 



U(R,4>) = U„\2\nR + g(^)], 

E(fl.0) = (5) 
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where Uq and Eq are constants. The functions and are 27r-periodic in because 
they are continuous. We normahze the surface density so that 

i^£/(0)d0=l. (6) 

This transfers the scahng of S to Eg (^-g-, GOl) and imphes Uq = ttGEq. Given /, Poisson's 
integral does not usually lead to a closed-form expression for g and one has to adopt a series 
solution. Let us assume that admits the Fourier expansion 

oo 

f{(j)) = 1 + J2 (Ci COS e(f) + Se sin £0) , (7) 
e=i 

where Q and S(, are constant coefficients. Using eqs (Al) and (A5) of ST96, we find the 
potentials corresponding to single Fourier modes and superpose the results. This gives 

g{(t>) = -2 J] Mcos^0+ -^sin^0j . (8) 



We prescribe / and determine g through the scheme mentioned above. Then, we attempt 
to solve the momentum equations for vr, v^, and the pressure p in terms of / and g. 

The velocity components can be derived from the stream function ^ as 

The scale-freeness of the model suggests the Ansatz 

*(i?,0)=Eoy^[lni? + P(0)], (10) 

where is a constant reference velocity that will become the circular velocity if the disk is 
axisymmetric (/ = 1, P = 0). We utilize eqs (5) through (10) and take the curl of eqs (2) 
and (3). This eliminates the pressure p and leaves us with the following ODE for P(0) 



p'" - y + j-,m'r - ff" + f]p' = pfifg' + 2/0, (11) 



where P — Uq/V^ is a dimensionless parameter and the prime denotes d/d(j). The solution 
P{4>) will result in the velocity field of our planar self-gravitating fiuids. Physical sohitions 
of (11) correspond to the condition v(P, 0) = v(P, + 27r), which ensures the continuity of 
velocity field. Hence, we seek for 27r-periodic solutions of P{4>). 
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3.2. Velocity Components 

Carrying out a change of dependent variable as (this is indeed Liouville's transformation 
explained in Appendix A) 

^' = ^(<^)y||, (12) 

eq. (11) reads 

w" + w^PmFo{(l>), Fo^fg' + 2f. (13) 
The particular solution of (13) is then obtained as 

i(;(0)=/3/(O) aosin0+6ocos0+y Fo(s) sin(0 - s)dsj , (14) 
where the (dimensionless) constant coefficients oq and bo are given by 

The function w{(f)), and therefore v, will be 27r-periodic if we impose the boundary conditions 
w{4> + 27t) = w{4>) and w'{4> + 27r) = w'(0) that lead to the necessary conditions 

j Fo(0)sin0d0= / Fo(0) cos0d0 = 0. (16) 
Jo Jo 

Using (12) and (9), the velocity components are readily determined as 
vr — — jiVrj) aosin0+6ocos0+ j Fo(s) sin(0 — s)dsj , 

from which one can obtain P{4>)- It is now straightforward to determine the pressure p from 
either (2) or (3) as 

^ ^ 1^ t'^^^'^^ " ^ ^^^^ 

where 

c^ = ;^[/9^(0)-l], (19) 

is the square of the sound speed (at the azimuth (p) and 

r f'^ 1 

^(0) — f{4>) {bo sin (f) — ao cos (/.) + 2/(0)- / Fo{s)cos{4>-s)ds\. (20) 
L JO 1 



-6- 



Sound waves are stable if > 0. i.e., p > 0. Define Amm — min{A(0) : < < 27r}. The 
necessary and sufficient condition for the positiveness of p will be 

< < ^min- (21) 

Apart from constraining the admissible values of this condition may put some restrictions 
on the model parameters. As it can be seen in (17), is maximum (minimum) where 
the surface density is minimum (maximum). This result has previously been used in the 
construction of DFs of stellar systems using Schwarzschild's orbit superposition method 
(e.g., K93; JZ). 

Our scale-free analytical solutions are general. They are independent of special angular 
variation of surface density and its corresponding potential. Our solutions can be used even 
in the presence of an external force that falls off as steep as R~^. Such a force can be that of 
a magnetic field or the gravitational influence of a dark matter halo. These effects can enter 
in our equations through Fo(0). 

The most important point related to our solutions is the existence of transonic flows. 

According to (19), it is possible to obtain entirely subsonic, transonic, and entirely supersonic 
flows depending on the choice of = 1/ -\/ PA^am in the range < 77 < 1. Let us define the 
Mach number 

= (22) 
which is a function of the azimuthal angle in our disks. For convenience, we use the variable 

MM = (23) 

for the classification of fiows. Disks with Mo{(p) > (A1o(0) < 0) for V0 G [0, 27r] are 
entirely supersonic (subsonic) and our steady-state models are credible. Flows for which 
A1o(0) switches sign along a streamline, are transonic. In such a circumstance, steady-state 
models become problematic as explained below. In principle, transition from subsonic to 
supersonic speeds (along a streamline) can take place smoothly, but the reverse is impossible 
without assuming a shock wave. Because fluid elements in the subsonic part of a streamline, 
have not information from the arrival of the particles in the supersonic upstream. Therefore, 
they suddenly find themselves collided with these particles. After a while, these interactions 
develop a shock wave and the speed of supersonic stream discontinuously decreases to a 
subsonic value. This process is inherently time-dependent and it increases entropy and 
generates heat. Thus, we must cautiously deal with any transonic solution. Galli et al. 
(2001) also argue that models with transonic fiows will most hkely be unstable. Despite these 
arguments, there is an exceptional case in nature where a supersonic fiowfield is isentropically 
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(and of course smoothly) compressed from supersonic to subsonic speeds. That is the famous 
conical flow detailed description of which can be found in Anderson (1990, chapter 10). In 
this very special case, again, a shock wave still presents but the sonic line does not coincide 
with the shock ray. The latter completely lies in the supersonic region. 

In any case, it seems reasonable to exclude transonic solutions from our results. A 
stability analysis will surely provide much insight into the problem. The instabilities of our 
non-axisymmetric disks will be resolved in a future paper. 

For axisymmetric disks, the parameter rj is obtained in terms of Mach number through 

In this case, M. is independent of 4> and the transition from subsonic to supersonic flows 
occurs at ?7 = l/-s/2- 



4. EXAMPLES 

We now apply our analytical tool to certain non-axisymmetric models for which we have 
either observational evidence or stellar dynamical counterparts. 



4.1. Simple Disks 



The simplest form of non-axisymmetric disks is obtained when a circular disk is per- 
turbed by a single Fourier mode as 

/(0) = l + /mcosm0, (25) 

where < |/^| < 1. The corresponding ^-dependence of the potential will become 

2 

g{(f))^ fmcosmcf). (26) 

m 

One can readily verify that the condition (16) is satisfied for all m > 1 and from (17) we find 

vr = PV^fm { , , sin 2m0 sin mcj)] . (27) 

V4m^ — 1 m + 1 / 

It is not difficult to show that A{(j)) is positive for our simple disks when < \fm\ < 1- These 
disks are analogous to the multiple- lobed disks of GOl. But they are not the same. The 
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difference can be traced back to the logic of solving the hydrodynamic equations. In fact, 
hydrodynamic equations for self-gravitating disks are four equations for the five unknowns 
E, U, Vr, and p. Thus, we are free to specify one of these variables, or put an extra 
constraint on them. Galli et al. (2001) follow the latter approach and assume an equation of 
state in the form p = a^S where a is the constant sound speed. They numerically solve their 
problem for S, v^, and vr, and obtain U through Poisson's equation. Our disks, however, 
are constructed relying on another approach. We prescribe the surface density E and find U 
using Poisson's integral. Then, we determine vr, v^, and p. Our disks, in contrast to GOl, 
are not isothermal [because the sound speed is not constant according to (19)]. Nonetheless, 
our solutions are exact and we never face the convergence problem as GOl. 

In Figure 1, we have surveyed the (77, fm) parameter space of our simple disks for ttG = 1, 
Eq = 1, and several choices of the wavenumber m. For each m, we obtain a curve (showed 
by solid, dashed, and dash-dotted lines) which is the boundary of models with transonic 
fiowfields. The region to the up side of each curve (indicated by the letter B) corresponds 
to the models that are transonic along a typical streamline. Below the boundary curve, we 
obtain entirely subsonic and entirely supersonic regions. Excluding transonic solutions (they 
are unphysical due to arguments in §3.2), highly non-axisymmetric disks can be constructed 
in the subsonic region as the wavenumber m is increased. The maximum elongation of the 
lopsided mode (m = 1) occurs for fm ~ 0.5 and r] ?a l/^/2. Whatever the value of m may be, 
the maximum value of 77 is 1/ \/2 for entirely subsonic models. Entirely supersonic models 
are approximately circular with < 0.2. 



There are two interesting families of bisymmetric disks for which self-consistent stellar 
dynamical models are known. They are the power-law and eUiptic disks (JZ). Disks whose 
equipotentials are similar concentric ellipses, are called the power-law disks. They are planar 
analogs of the logarithmic eUipsoidal models (e.g., Binney 1981; Evans, CaroUo & de Zeeuw 
2000). The power- law disks with fiat rotation curves take the potential of the form (5) with 



where q is the axial ratio of equipotentials. The angular dependence of the surface density 
corresponding to (28) has been derived in JZ as the following convergent series 



4.2. Bisymmetric Disks 



5f(0) = ln(g^ cos^ (p + sin^ 0), 



(28) 



00 




(29) 



even ^>0 
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(1 - >i)"s. 



kl 



)n-l 



n=l 



n-1 
m=0 



2n 



with A; = 2n, i = 
have defined (3m '' 



n + m, and j = n — m. In this relation, is the Kronecker delta and we 
= if n — m is odd, and 



g(n) ^ 



n 

n — m\ 



m ^ 0, 



(30) 



n 



m 



0, 



if n — m is even. 

When the isodensity contours are similar concentric ellipses, we attain elliptic disks with 

/(0) = (g^cos2 + sin2 0)-V2, (31) 

where q is the axis ratio of isodensity contours. Stellar elliptic disks were constructed (for 
the first time) in K93 and then were examined for self-consistency by the curvature condition 
in JZ. The potential function associated with (31) was determined by Evans & de Zeeuw 
(1992, eq. 5.2) in terms of an indefinite integral, but a practical way is to generate it using 
the scheme presented in §3.1. We find 



oo p 

U{R,(t)) = 27rG'Eo(lni?- y cos 

even €>0 

(1 - qYV{\ + k)V{\ + n - k)(3^^\5u + M 



(32) 



n=l fc=0 m=0 



27rA;!(n- A;)! 



with i = 2k — n + m and j = 2k — n — m. All the series introduced so far, rapidly converge 
using only the first few terms. 

For instance, we have computed the isodensity contours and streamlines of our bisym- 
metric disks for q = 0.5, nG = 1, and Sq = 1, which result in Amin = 0.854 and Amin = 3.041 
for the power-law and elliptic disks, respectively. Figure 2 shows the results of our calcu- 
lations for several choices of V^. Solid lines are the isodensity contours and dashed fines 
are streamlines. For the small values of V^, the pressure increases and dominates centrifu- 
gal forces. Therefore, streamfines become highly dimpled in both disks (Panels 2a and 2c). 
As the value of is increased (V^ — > UoA^iri, i-e., (3 — > l/ylmin), centrifugal forces over- 
whelm the pressure forces and keep the fiuid particles away from the center. This causes 
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streamlines to become rounder (Panels 2b and 2d). Numerical computations show that ^(0) 
(0 < < 2n) is positive for the power-law and elliptic disks in the range < g < 1. There- 
fore, it is always possible to find a physical value for V^. Figure 3 shows the (77, q) parameter 
space of the power-law and elliptic disks. The region filled by dots corresponds to models 
with transonic fiowfields. They are most likely unphysical models as we discussed in §3.2. 
Hence, we exclude them from our solution space. As Figure 3 shows, supersonic elliptic disks 
can be more elongated than supersonic power-law disks. In both cases the maximum value 
of 77 for subsonic flows is equal to 1/V2. This is a generic property of our self-gravitating 
gaseous disks inherited from simple disks. 

Another remarkable point is the limit of the axis ratio q for plausible models (transonic 
models excluded). The lower bound of q is almost 0.4 in eUiptic disks but in the case of 
power-law disks one can see a narrow subsonic region with q < 0.08. In the limit of g — > 0, 
these highly elongated models converge to a needle. The results for eUiptic disks are similar 
to the findings of JZ for stellar disks. In fact, stellar elliptic disks with q ^ 0.5 are non- 
self-consistent. On the evidence of Figure 3a, gaseous power-law disks exhibit significantly 
different features compared to their stellar counterparts that cannot be self-consistent if 
q < 0.707 (JZ). 

Our calculations show that the gas pressure of elliptic disks approximately takes an 
axisymmetric distribution when is small. As a result, the line-of-sight velocity dispersion 
a = (wfos) — (wios)^ will inversely be proportional to /(0). This property does not depend on 
the choice of q and it predicts a correlation between the photometric and spectroscopic data 
of elliptic disks when centrifugal forces are dominated by internal stresses. 

5. IMPLICATIONS FOR THE KINEMATICS OF THE LMC 

In his recent study of the LMC, van der Marel (2001) has shown that the LMC disk is 
elongated and its isophotal lines are well fitted by ellipses. Since the LMC has an approxi- 
mately flat rotation curve, our elliptic disks (presented in §4.2) can be applied to model the 
LMC disk. 

The most recent studies of the kinematics of gas and discrete tracers (like carbon stars) 
of the LMC (Kim et al. 1998; Alves & Nelson 2000) show that the line of maximum velocity 
gradient is in the range 143° < Omax < 183°, which differs from the angular position of the 
line of nodes (6 = 122°.5) by 20° < Q^^ - 9 < 60°. Noting the fact that the intrinsic 
ellipticity of the LMC disk is e = 0.312 in outer regions, van der Marel (2001) argues that 
A© = ©max — © ~ 20° is possible but it is hard to believe that A© can be as large as 
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60°. He ascribes the excess of computed misalignment to the center-of-mass and rigid-body 

motions of the LMC. A more accurate kinematic analysis of the LMC predicts a value of 
A0 as small as 7.4° (van der Marel 2002, private communication). To provide a theoretical 
framework for these observations, we have modeled the outer region of the LMC using our 
elliptic disks. 

We have calculated the position angle of the maximum line-of-sight velocity of our 
elliptic disks, given the viewing angles © = 122°. 5 and i = 34°. 7 that were determined 
by van der Marel & Cioni (2001). Our results have been shown in Figure 2. This figure 
illustrates the variation of A© versus r] for several choices of e = 1 — q including e = 0.312, 
which is the intrinsic ellipticity of the LMC at large radii. A discontinuity occurs in the graph 
of A© (let us say at r] = because of a sudden change in the morphology of streamlines. 
For 77 < rjcr, streamlines are dimpled. The graph for e = 0.312 (filled squares) shows that 
A© varies between 0° and 30° when streamlines are not dimpled. Interestingly, this result is 
in harmony with the predictions of van der Marel (2001). It is remarked that for e = 0.312 
we obtain = 0.36. As 77 ^ 1, streamlines become almost circular independent of the 
value of e. In such a circumstance, the position angle of the maximum line-of-sight velocity 
coincides with that of the line of nodes, i.e. A© 0, as expected. Based on the available 
observational data and our results for e = 0.312, we conclude that our scale-free elliptic disks 
(with flat rotation curves) can explain the kinematics of the outer region of the LMC if 77 
is chosen in the range 0.36 < rj < 0.72. The upper bound on rj comes from the fact that 
our transonic solutions are not reliable answers of steady-state equations. For A© = 7.4°, 
we find 77 0.54, which is well within the acceptable range. Moreover, models with 77 ^ 1 
(they are supersonic) lead to very high fluctuations of the pressure, and therefore, velocity 
dispersion. This does not seem to be relevant to a stable model. We also rule out the 
possibility of 77 < rja-, for it gives A© < 0, which is not observed. Furthermore, the condition 
p — > at the boundary of the LMC disk is inconsistent with a pressure-dominated model. 



6. CONCLUSIONS 

Exact equilibrium solutions of self-gravitating systems are rare. Two classes of such 
solutions were found by Medvedev & Narayan (2000) for three-dimensional, axisymmet- 
ric isothermal systems. A generalization of these solutions for non-axisymmetric rotating 
systems was then introduced by Shadmehri & Ghanbari (2001). Although the density distri- 
bution of their models is non-axisymmetric, the velocity field has only the toroidal compo- 
nent. This makes their models too simple to be used in the study of fully three-dimensional 
non-axisymmetric distributions of matter. 



- 12 - 



In a systematic search for non-axisymmetric equilibrium solutions, we started our study 
from gaseous disks because they inherit many properties of three-dimensional systems. Using 
the benefit of the existence of a stream function, we were successful in reducing the governing 
equations of scale- free models to an incomplete linear ODE for the angular part of the stream 
function. It then became clear that Liouville's transformation can reduce the obtained ODE 
to a second-order inhomogeneous ODE, which is similar to that of a forced harmonic oscilla- 
tor. This was striking because it led us to derive closed-form solutions for the radial velocity 
component, stream function and the pressure p. In fluid mechanics, transport equations are 
written for mass, momentum and energy. While the pressure occurs in governing equations, 
there is no extra independent equation in order to determine p. The problem is rather diffi- 
cult in the presence of an anisotropic stress tensor, which is encountered in turbulence and 
Jeans' equations of stellar hydrodynamics. Taking the curl of momentum equations is a trick 
that temporarily removes the pressure from our calculations of the stream function. But, one 
should take care of the sign of p once it is eventually extracted from momentum equations. 
Only a positive pressure distribution is physical because a fluid medium can resist against a 
compressive stress not a tensional one. 

Non-axisymmetric stellar disks have necessarily anisotropic DFs and pressure tensors. 
This point is the main difference between the physics of a steady-state gaseous disk and its 
stellar counterpart described by Jeans' equations. Nevertheless, isotropic pressure distribu- 
tion is not an obstacle for the existence of non-axisymmetric, self-gravitating gaseous disks 
as we showed in §4. 

Our scale-free elliptic disks were used to model the outer part of the LMC disk. The main 
restriction of our modehng is the isotropy of the stress tensor. We speculate that our results 
will not be changed in a significant way if the model is stellar with an anisotropic pressure. 
The other important issue is the perturbation of the Galaxy. The tidal field of the Galaxy 
induces a force whose i?-distribution is totally different from that of a logarithmic potential 
field. However, if the tidal field of the Galaxy was dominant, the rotation curve of the LMC 
would show substantial deviation from flatness. But, we know that the rotation curve of the 
LMC is almost flat. Therefore, our isolated models provide acceptable zeroth- order so\\iiions. 
Such a modeling of the LMC is a quick route for interpreting the observational data. It is 
obvious that taking the tidal effects into account, will destroy the enjoyable scale-free nature 
of our disks. 

We thank Tim de Zeeuw for useful comments and drawing our attention to the pos- 
sible application of our findings to the LMC. We are indebted to an anonymous referee 
for enlightening suggestions that helped us to substantially improve the presentation of our 
results. 
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LIOUVILLE'S TRANSFORMATION 



Consider the following non-autonomous, linear, ordinary differential equation for Z in 
terms of t: 

d^Z „ _ dZ 



+ Qi(t)— + g2(t)^ = goW. (Ai) 



di2 ' dt 

Using the transformation (Bellman 1966, page 82; Verhulst 1996, problem 8-2) 



Z{t) = w{t) exp (^-^ ' 



which has been introduced by Liouville, eq. (Al) becomes 



+ [Q2 - - = QoWexp {llQiir)dT^ ■ (A3) 



di2 

The study of this equation is easier than (Al). 

We assume P' as the dependent variable. This puts eq. (11) in the form of (Al). We 
then apply eq. (A2) and obtain eqs (12) and (13). 
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Fig. 1. — The {r],fm) parameter space of simple disks. Solid, dashed, and dash-dotted 
lines respectively show the boundary of transonic models for m = 1, m = 2, and m = 3. 
More elongated subsonic models are supported as the wavenumber m is increased. Entirely 
supersonic models can exist for a limited ranges of parameters. Such models are almost 
axisymmetric. 

Fig. 2. — Isodensity contours (sohd hnes) and streamhnes (dashed lines) of the power- law 
and elliptic disks with fiat rotation curves for Eq = 1, ttG = 1, and q = 0.5. (a) Power-law 
disk with = y/A~/2 = 0.462. (6) Power-law disk with = A^/A~/5 = 0.739. (c) 
Elhptic disk with - = 0.581. (d) Elliptic disk with V4, = 4:^/A~/5 = 1.395. 

Fig. 3. — The (77, q) parameter space of the power-law (Panel a) and elliptic disks (Panel 
b). Dotted region corresponds to models with transonic fiowfields, which are most likely 
unphysical. Elliptic disks allow for more elongated supersonic models than the power-law 
disks. The maximum value of r] taken by subsonic models is equal to 1/a/2 in both disks. 

Fig. 4. — Variation of A© = ©max — © versus 77 = i/VPA^ for elliptic disks, given the 
viewing angles © = 122°. 5 and i — 34°. 7. Circles, filled squares, and deltas correspond to 
the eUipticities e = 0.15, 0.312, and 0.6, respectively. The existing discontinuities in the data 
sets are due to a sudden change in the morphology of streamlines. This happens at a critical 
T) — rjcr- For T) < TjcT, streamlines are dimpled. 
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